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ABSTRACT 

We determine the gravitational response to a diffuse source, in a locally 
de Sitter background, of a class of theories which modify the Einstein-Hilbert 
action by adding a term proportional to an inverse power of the Ricci scalar. 
We find a linearly growing force which is not phenomenologically acceptable. 
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1 Introduction 



Observations of Type la Supernovae give compelling evidence that the uni- 
verse is entering a phase of acceleration pQ. This means the current energy 
profile is dominated by an unknown source contributing negative pressure 
which has become known as "dark energy". When various recent data sets 
are combined, the fraction of dark energy present in the universe is deter- 
mined to have the value, Q\ ~ 0.73 [2 . The number and quality of the most 
recent measurements [T] leave little doubt that we are facing a real effect 
which must be explained. 

There has been no lack of theoretical effort to account for and describe 
dark energy origins and dynamics None of the suggestions is without 

problems. A bare cosmological constant will work, but one has to understand 
both why it is more than 120 orders of magnitude smaller than its seemingly 
natural scale, and why it has just achieved dominance in the current epoch 
[miHI]. Quintessence based on a scalar field will also work [12] -[El but one 
must understand why it is homogeneous [17] and again why it has achieved 
dominance now. Long range forces ^H] and even quantum effects ^H] have 
also been suggested. 

A recent paper by Carroll, Duvvuri, Trodden, and Turner proposed a 
purely gravitational approach [20] ■ Late time acceleration is achieved by 
considering a subset of nonlinear gravity theories in which a function of the 
Ricci scalar is added to the usual Einstein-Hilbert action, 



where p > and \i is an a priori unknown parameter. It has been argued that 
such inverse powers of R may be closely related with braneworlds and string 
theory [21]. By considering the standard Friedmann- Robertson- Walker cos- 
mology, Carroll, Duvvuri, Trodden and Turner showed that if /i ~ 10 _33 eV 
(the inverse age of the universe) the equation of state parameter could fall 
within the range —1 < w c s < —2/3 at late times. 

Although higher derivatives typically bring negative energy degrees of 
freedom, endowing the Lagrangian with nonlinear functions of the Ricci 




(1) 




(2) 



1 



scalar can sometimes be acceptable (221- This will only give rise to a sin- 
gle, spin zero higher derivative degree of freedom. Since the lower deriva- 
tive spin zero field is a constrained, negative energy degree of freedom (the 
Newtonian potential), its higher derivative counterpart can sometimes carry 
positive energy 

Recent papers have considered different aspects of the model. Dick con- 
sidered the Newtonian limit in perturbation theory about a maximally sym- 
metric background [22] • Dolgov and Kawasaki discovered an apparently fatal 
instability in the interior of a matter distribution [21]. However, Nojiri and 
Odintsov have shown that an R 2 can be added to the Lagrangian without 
changing the cosmological solution, and that the coefficient of this term can 
be chosen to enormously increase the time constant of the interior instability 
|25j . Meng and Wang have explored perturbative corrections to cosmology 
[26J. Others have drawn connections with a special class of scalar-tensor 
theories |27|-|2B|. 

Our own work concerns the response to a diffuse spherical matter source 
after the epoch of acceleration has set in. The procedure will be to solve 
for the perturbed Ricci scalar, whence we determine the gravitational force 
carried by the trace of the metric perturbation. We constrain the matter 
distribution to have the property that its rate of gravitational collapse is 
identical to the rate of spacetime expansion, thereby fixing the physical radius 
of the distribution to a constant value. Further, we impose the condition that 
inside the matter distribution the density is low enough to justifiably employ 
a locally de Sitter background, in which case the Ricci scalar can be solved 
exactly and remains constant. 

Our de Sitter background was considered by Carroll et al. [20 a and was 
found to be unstable. However, the decay time ~ /x -1 is far too long to create 
any practical concern. The next section describes the calculation and shows 
the solution outside the matter distribution to grow linearly. We conclude 
by remarking on the implications of this result. 

2 The Gravitational Response 

We shall consider a gravitational action parameterized by p > 0, 

S[g) = ^ J d A xV^ [R - fi 2{p+1) R~ p ] . (3) 
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(We employ a spacelike metric with Ricci tensor R pv = r p — ^ P PP)U + 
r p r°j, — T P V(T Y C .) Functionally varying with respect to the metric and 
setting it equal to the matter stress energy tensor leads to the equations of 
motion, 

[1 +p^R-^} R pv -\\\- ^R-^] Rg, u 

+ Pfi 2(p+1) (g^D - D^D V )R-^ = SttGT^. (4) 

is the covariant derivative and □ = (—g)~ 1 ^ 2 d ll (y/—gg^ u d v ) is the covari- 
ant d'Alembertian. 

Although one must really solve all components of the field equations (@J) 
we can get an important part of the gravitational response by simply taking 
the trace. We shall also restrict to p — 1 for simplicity. Inside the matter 
distribution the trace equation is, 

-R + 3^ + 3/i 4 D-^ = SnGg^T^ = T. (5) 

(Note that T is negative.) Normally, one would expect the matter stress en- 
ergy to be redshifted by powers of the scale factor in an expanding universe. 
However, recall that this matter distribution possesses a rate of gravitational 
collapse equal to the rate of universal expansion, and thus T remains con- 
stant. Since our matter source is also diffuse, we may perturb around a 
locally de Sitter space. For the interior solution, we are able to solve for R 
exactly using equation (|SJ) for the case T is constant and D p R = 0, 



Rm — 
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Obtaining de Sitter background obviously selects the negative root. Further, 
we concentrate on the situation |T| -C fi 2 , 

R in = vV - y + • • " ■ (7) 

Outside the matter source we perturb around the de Sitter vacuum solu- 
tion, 

Ro»t = \/3/i 2 + SR. (8) 
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Substituting (jSJ) into (0) and expanding to first order in 5R yields the equa- 
tion defining the Ricci scalar correction, 

D5R(x) + V3/j 2 5R(x) = 0. (9) 

In our locally de Sitter background the invariant length element is, 

ds 2 = -dt 2 + a 2 (t)dx ■ dx, (10) 

with a(t) having the property, 

d 

H = — — constant. (11) 
a 

We can relate the Hubble constant H to the parameter fi via the vacuum 
Ricci scalar, 

R = 12H 2 = V3/J 2 . (12) 
Identifying □ = a~ 3 d p (a 3 g p(r d (T ), we expand ©, 

[d 2 - 3Hd + 12H 2 ] 6R(t, x) = 0, (13) 

where d 2 = — d 2 + a _2 V 2 . It is evident from (|T3*j) that the frequency term 
has the wrong sign for stability [20] ■ However, since the decay time is pro- 
portional to 1/H, we may safely ignore this issue. 

Seeking a solution of the form 5R = 5R(Ha\\x\\) allows us to convert (|13|) 
into an ordinary differential equation, 

5R = 0, (14) 
where y = Ha\\x\\. To solve this equation we try a series of the form, 

oo 
n=0 

Substituting this series into (|T3*|) yields a solution with a = 0, 

'•'"'^^RniTfrw (1B) 



d 



d 



dy 



dy 



4 



and a solution with a = 



-1, 



_57^ 

4 - 



_57^ 
4 ) 



(2y) 



2 ci 



#)r(i + ^f) (2^)! 
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Both solutions converge for < y < 1. Both also have a logarithmic singu- 
larity at y — 1, which corresponds to the Hubble radius. We can therefore 
employ them quite reliably within the visible universe. 
The solution we seek is a linear combination, 



<Ji2(y)=/3i/o(y)+ft/-i(y), 



(18) 



whose coefficients are determined by the requirements that SR(y) and its 
first derivative are continuous at the boundary of the matter distribution. 
We employ a spherically symmetric distribution of matter, centered on the 
co-moving origin. If the matter distribution collapses at the same rate as the 
expansion of the universe, its physical radius is a constant we call p. (This 
means that the co-moving coordinate radius is p/a(t).) If the total mass of 
the distribution is M we can identify T as the constant, 



r 



~4 



6GM 



(19) 



In terms of our variable y = Ha(t)\\x\\, the boundary of the matter distri- 
bution is at uq = Hp. Demanding continuity of the Ricci scalar and its first 
derivative at Uq gives the following result for the combination coefficients of 
the exterior solution (JTHJ), 
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where a prime represents the derivative with respect to the argument. 

We are now in a position to calculate the gravitational force carried by the 
trace of the graviton field. The metric perturbation modifies the invariant 
length element as follows, 



ds 2 



[1 - h 00 )dt 2 + 2a(t)h 0i dtdx i + a 2 (t)(5ij + h ij )dx i dx j . 



(22) 
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Further defining h = —h 00 + ha and imposing the gauge condition, 



allows us to express the Ricci scalar in terms of h, 



5R - A - ( -(V 2 // -\ IHOJi 



(23) 



(24) 



(Recall that we define = (<9 t ,a X V).) Assuming ft, = h(y) as we did for 
5R gives the equation for the gravitational force carried by h, 



(l/ 2 -l)^ + -(5y 2 -2) 
dy y 



h'(y) 



25R(y) 
H 2 



(25) 



The solution to (|25|) is 

/»'(*) = 



y 2 (l - y 2 f' 2 J 



dy'y' 2 (l- y y/ 



1/2 «W 



# 2 



(26) 



At this point it is useful to consider the y values which are relevant. The 
Hubble radius corresponds to y = 1, whereas the typical distance between 
galaxies corresponds to about y = 1CT 4 , and a typical galaxy radius would 
be about y = 1CT 6 . We are therefore quite justified in assuming that y <C 1, 
and in specializing to the case of y <C y <C 1. Now consider the series 
expansions, 



f (y) = l-2y 2 + \y* + O( y 6 
5 



1 - 7y 2 + ^ + 0(y e ) 



3MG 2 , 
\2MGyl 



0{yl 



(27) 
(28) 



We see first that |/3 2 | <C /3i — which means SR(y) ~ Pifo(y) — and second, 
that /o(y) ~ 1 - - which implies $R{y) ~ — T/2. This means that the 
integrand in (J2H)) fails to fall off for y > y , so the integral continues to grow 
outside the boundary of the matter distribution. For small ?/ > j/ we have, 



h'(y) 



2GM 



y + o( y 3 ). 



(29) 
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To see that this linear growth is a phenomenological disaster it suffices 
to compare (|29|) with the result that would follow for the same matter dis- 
tribution, in the same locally de Sitter background, if the theory of gravity 
had been general relativity with a positive cosmological constant A = 3H 2 . 
In that case SR(y) = —T9(y — y) and, for y > y , the integral in (j2Tlj) gives, 

h '(v) nB = AZJ2 Z ^{axcsinfoo) - y (l - 2yg) J I - y 2 \. (30) 

gr 4H z y z {l — y z )2 < v > 

= J^l + OW . (31) 

y 

The linear force law (|29j) of modified gravity is stronger by a factor of |(^) 3 - 
For the force between two galaxies this factor would be about a million. 



3 Conclusion and Remarks 

We have determined the gravitational response to a diffuse matter source 
in a locally de Sitter background. Our result is the leading order result 
in the expansion variable y, the fractional Hubble distance. Equation (|29|) 
clearly forces us to disregard the class of theories considered here (J2J) when 
compared to GR with a cosmological constant (for example, the correction 
to the gravitational force between the Milky Way and Andromeda increases 
by six orders of magnitude). 

The two assumptions made in our analysis were: 

• the matter distribution is gravitationally bound, 

• the matter distribution has a mean stress energy |T| < /i 2 . 

The second of these assumptions can be viewed rather flexibly if interested 
only in phenomenological implications. Regardless of whether it is satisfied, 
we still would expect a linearly growing response far from the source. To 
see this, recall that the dominant piece of the solution, fo(y), from equation 
()18|) remains constant and approximately equal to one for many orders of 
magnitude (for instance, /o(10~ 8 ) — /o(10~ 3 ) ~ 1CT 6 ). Therefore, although 
the exterior solution would not be very reliable near the matter source, we 
can be confident that at cosmic or even intergalactic scales perturbing about 
de Sitter becomes appropriate and a growing solution would still be observed. 
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This analysis was performed for p = 1, but of course nothing restricts us 
from considering arbitrary powers of the inverse Ricci scalar. To no surprise, 
however, varying the power only changes the coefficient of the gravitational 
force leaving its qualitative behavior alone. The instability found by Dolgov 
and Kawasaki [24J and the growing solution calculated in this work seem 
to preclude all such theories phenomeno logically. The two problems seem 
to complement one another because either problem could be avoided by the 
addition of an R 2 term, which would not alter the cosmological solution 
|25j . However, avoiding the interior instability seems to require the R 2 term 
to have a large coefficient, whereas avoiding the exterior growth requires a 
smaller value [2*U] . 

None of these issues diminishes the importance that should be placed 
on considering novel approaches to understanding the dark energy problem. 
It is the responsibility of both theorists and experimentalists to construct 
and constrain candidate theories, and it is truly an exciting epoch of human 
investigation for which we are just beginning to acquire these capabilities. 
Greater freedom can be obtained by adding different powers of R. (Note that 
this generally alters the cosmological solution.) Although such models seem 
epicyclic when considered as modifications of gravity, the same would not be 
true if they were to arise from fundamental theory. For example, it can be 
shown that the brane-world scenario of Dvali, Gabadadze and Shifman [3*U] 
avoids both the interior instability and the linearly growing force law [3Tj . 
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